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Abstract. The problem of determination of the orbital velocity of an astrometric satellite from 
its own observational data is studied. It is well known that data processing of microarcsecond- 
level astrometric observations imposes very stringent requirements on the accuracy of the orbital 
velocity of the satellite (a velocity correction of 1.45 mm/s implies an aberrational correction of 
1 ^tas). Because of a number of degeneracies the orbital velocity cannot be fully restored from 
observations provided by the satellite. Seven constraints that must be applied on a velocity pa- 
rameterization are discussed and formulated mathematically. It is shown what part of velocity 
can be recovered from astrometric data by a combined fit of both velocity parameters and astro- 
metric parameters of the sources. Numerical simulations show that, with the seven constraints 
applied, the velocity and astrometric parameters can be reliably estimated from observational 
data. It is also argued that the idea to improve the velocity of an astrometric satellite from its 
own observational data is only useful if the a priori information on the orbital velocity justifies 
the applicability of the velocity constraints. The proposed model takes into account only transla- 
tional motion of the satellite and ignores any satellite-specific parameters. Therefore, the results 
of this study are equally applicable to both scanning missions similar to Gaia, and pointing ones 
like SIM, provided that enough sources were observed sufficiently uniformly. 
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1. Introduction 

The astrometric accuracy of several microarcseconds announced for three space mis- 
sions (Gaia, SIM and Jasmine) implies that the velocity of the satellite should be known 
within several mm/s. Such a high precision is a challenge to standard orbit determination 
techniques. 

Another way to proceed was discussed in the Gaia community since 2001 and formu- 



lated in a written form by Klioner (2005) The idea is to use Gaia's own astrometric 



data to fit a correction to the Gaia velocity. The velocity correction 5\ is the difference 
between the real velocity of the satellite and its ephemeris velocity available a priori. 
A straightforward theoretical analysis of this idea (Butkevich 2006) demonstrated that 
velocity can be determined from observations with the required precision, provided that 
all other parameters are exactly known. Some criticism against this approach has been 



formulated by Bastian (2004a) who argued that the velocity obtained from observation 
will strongly correlate with the source parameters. A detailed theoretical exposition of 
this problem has been done by |Klioner fc Butkevich (2007) who explicitly demonstrated 



that the velocity correction to be fitted from the data must satisfy some constraints, 
t On leave from Pulkovo Observatory, 196140 Saint-Petersburg, Russia. 
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2. What part of observer's velocity can be restored from astrometric 
data? 

Let us consider the case when the velocity correction has a constant component: 

5v = vo = const , (2.1) 

then this additional velocity leads to the following aberrational correction 

5u « -u x [u x vo/c] . (2.2) 

This correction depends on the stellar position u, but for each star it is constant. Such 
a signal in <5v is equivalent to a constant change of positions, which cannot be detected 
from observations since it cannot be distinguished from their real change. 
Similarly, for a velocity correction changing linearly with time, 

5v = a t a = const , (2.3) 

we have 

5u(t)«/V, (2.4) 

where 

(Ji = — u x [u x a /c] . (2.5) 

This correction is equivalent to an additional constant proper motion for each star. Such 
a correction again cannot be detected from observations since it cannot be distinguished 
from a real change of the proper motion parameters for each source. 

If <5v is exactly proportional to the barycentric position of the satellite r (t): 

5v (t) = aor (t) a = const , (2.6) 

the corresponding first-order aberrational correction reads 

8u w — u x [u x a r/c] . (2.7) 

On the other hand, the aberrational effect caused by a global offset of parallaxes Sn is 

Su = u x [u x Snr/AU] . (2.8) 

These effects are indistinguishable provided that 

Sir = -a AU/c. (2.9) 

Thus the problem has seven free parameters («o and six components of Vo and ao), 
which correlate with some astrometric information, i. e. it has seven degrees of freedom. 
This rank deficiency makes the direct determination of velocity impossible. It can be, 
however, demonstrated that the degeneracy can be eliminated if the following constraints 
would be imposed onto the solution: 



i 

I 



T 

5v (t) dt = to remove Sv = v , (2.10) 



o 

T 



/ 



(t-T/2)5v(t)dt = to remove 6v = a t , (2.11) 

^ v (t) r (t) dt = Q t0 rcmovc ,5 V = a r (2.12) 
r(t)| 2 
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Determination of the barycentric velocity from observations 
Table 1. Singular values of the normal matrix. 
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n 1:7 8 : max 

cr n 9.6 • 1(T 8 : 9.0 ■ 1(T 6 2.1 : 1045 



Table 2. Error in parameters. 



Parameter 


Oi S [la. 


[15 7T 


Sv x Sv y 8Vz 


Error 


1.6 /ias 1.5 /^as 1.3 /ias/yr 


1.2 fias/yr 1.5 /ias 1.1 


3 mm/s 1.5 mm/s 1.6 mm/s 



These constraints guarantee that the solution does not contain the relevant signals in 
the sense of least-squares. 

Although it can hardly be proved analytically that the problem has no other degrees 
of freedom, it can be checked numerically. We calculated singular value decomposition 
(SVD) of a normal matrix and found only seven small singular values shown in Table [TJ 
This fact states that no other degrees of freedom exist. 



3. Legitimacy of the constrained velocity 

The constraints may be only applied when a priori accuracy of the ephemeris guar- 
antees that no signal of given kinds can exist in real velocity, or, strictly speaking the 
signal is so small that any effect due to it can be completely neglected. 

The real 5v roa i velocity correction can be represented as a sum of the two different 
components 

<5v rca i = <5v + R (t) , (3.1) 

where 5v is the component that can be fitted and R (t) is the component violating the 
seven constraints. The fitted component (5v is useful if and only if the uncertainty of the 
ephemeris velocity is such that it guarantees that 

|R(t)|<e (3.2) 

at any instant of time. Here e is a required velocity accuracy (for Gaia, for example, 
e = 1 mm/s). 

Fortunately, this can be demonstrated for Gaia (Klioner & Butkevich, 2007) but cannot 
be guaranteed for other missions. 



4. Results of the numerical simulations 

To study the problem of the velocity determination numerically, we have implemented 
a simple simulator of Gaia observations - Dresden Gaia Simulator (DGS) - that includes 
many of the basic features of the real mission. Table [2] shows the results of a simula- 
tion run with 2048 stars covering 5 years of observations (Butkevich & Klio ner 2007ap . 
The accuracy of an individual observation was chosen to be 30 /Ltas. The small errors 
found suggest that the constrained solution allows one to achieve a precise and reliable 
determination of velocity and source parameters. 
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Figure 1. Dependence of the velocity accuracy (the descending curve) and the number of 
stars in the astrometric solution (the ascending curve) on the limiting magnitude. 

5. Assessment of the accuracy of velocity determination 

Our simulation has one serious drawback - it can only handle limited datasets and 
never approaches the data volume close to the expected Gaia parameters. A straightfor- 
ward estimation can be made using simple statistical considerations (Butkevich & Klioncr 2 007b [) . 
One of the critical parameters in the velocity determination is the number of stars used 
in the solution, which depends on the apparent magnitude cut-off. The estimated ve- 
locity accuracy together with relevant star counts are shown in Fig. [T] for the limiting 
V magnitudes. Besides the number of stars, the accuracy also depends on the temporal 
resolution of the fitted velocity correction. The accuracy obviously degrades when a finer 
resolution is used. The time scale of velocity variations was chosen to be 6 hours, close to 
the rotation period of Gaia. We may conclude from the obtained estimates that at least 
10 6 stars are needed to obtain velocity with an accuracy of 1 mm/s. 

6. Satellite specific parameters 

The discussed model takes into account only translational motion of the satellite and 
ignores any satellite-specific calibration parameters (e. g. attitude parameters). The sit- 
uation may become more complicated when those other parameters are also considered. 
Our analysis shows that for Gaia scientifically important parameters can be successfully 
restored even in this case. 
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